Delay-Time and Thermopower Distributions 
at the Spectrum Edges of a Chaotic Scatterer 
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We study chaotic scattering outside tiie wide band limit, as tlie Fermi energy Bp approaches 
the band edges Eb of a one-dimensional lattice embedding a scattering region of M sites. The 
Hamiltonian Hm of this region is taken from the Cauchy orthogonal ensemble. The scattering is 
chaotic at Ep if the average level density per site of Hm at Ep describes a semi-circle as Ep varies 
inside the conduction band. The edges of this semi-circle coincide with the band edges Eb- We 
show that the delay-time and thermopower distributions differ near the edges from the universal 
expressions valid in the bulk. To obtain the asymptotic universal forms of these edge distributions, 
one must keep constant the energy distance Ep — Eb measured in unit of the same energy scale 
oc M~^^^ which is used for rescaling the energy level spacings at the spectrum edges of large 
Gaussian matrices. In particular the delay-time and the thermopower have the same universal edge 
distributions for arbitrary M as those for an M = 2 scatterer, which we obtain analytically. 

PACS numbers: 73.23.-b 42.25.Bs, 72.20.Pa 05.45.Mt 



Since nano- engineering makes it possible to fabricate 
devices which can be used for harvesting energy (Seebeck 
effect) or for cooling (Peltier effect) at the nanoscales, 
to study thermoelectric conversion beyond the semi- 
classical Boltzmann limit has become an important chal- 
lenge. In the mesoscopic limit, quantum interferences 
can induce large fluctuations and the thermoelectric con- 
version can be mainly due to a large fluctuation around 
a zero average. This requires the knowledge of the whole 
distributions instead of simple averages. Such mesoscopic 
fluctuations of the thermopower have been observed [ij 
in quantum dots as one varies either the shape of the dot 
with an electrostatic gate or the magnetic flux thread- 
ing the dot. Scattering theory, combined with random 
matrix theory (RMT), allows one to obtain Q these dis- 
tributions as long as the device behaves as an elastic scat- 
terer (the low temperature limit). 

Theoretical evaluation of the distribution of the ther- 
mopower Sk at a low temperature T is complicated 
by the fact that one needs to know the distribution of 
not only the scattering matrix S but also its energy 
derivative at the Fermi energy Ep, which can be ob- 
tained from the Wigner-Smith time-delay matrix 0-0] 
Q = —ihS~^dS/dE. The universal distributions of Q, 
and hence of Sk, have been obtained previously 0, but 
only in the wide band limit (WBL), restricted to the hulk 
of a wide conduction band. In that case, this is the bulk 
of the spectrum of the used chaotic scatterer which is 
probed at Ep. In the microscopic model which we study 



exactly, one can consider not only the bulk of the con- 
duction band, but also its edges where this is not bulk 
of the spectrum, but its edges which are probed at Ep. 
This is particularly interesting, since the universality of 
a distribution near the spectrum-edges can be quite dif- 
ferent from that in the bulk, as is well-known from the 
Tracy- Widom vs Wigner nearest-neighbor spacing distri- 
butions [^-Q. In this letter we show that the distribu- 
tions of Q and Sk near the edges are indeed different, and 
that the edge distributions give rise to a new asymptotic 
universality when the Tracy- Widom scaling is adopted. 
Furthermore, the results concerning Q which we obtain 
near the band edges could be relevant for other waves Q 
(electromagnetic, acoustic, . . .) than electron waves. 

Let us consider a chaotic cavity opened to leads via 
single mode quantum point contacts. We assume time 
reversal and spin rotation symmetry, such that the scat- 
tering matrix 5* is a 2 x 2 unitary symmetric matrix, di- 
agonalizable by an orthogonal transformation. The gen- 
eralization to other symmetries (unitary and symplectic) 
is straightforward, and will not be discussed here. S{Ep) 
can be assumed to be totally random at the Fermi en- 
ergy Ep, the probability P{S) of finding S inside an in- 
finitesimal volume element dS being constant. In other 
words, S(Ep) belongs to the Circular Orthogonal Ensem- 
ble (COE) \l(f\. Such an assumption can be justified if 
the corresponding classical trajectories are chaotic (9llT]|. 
The two eigenvalues e**^ of S{Ep) exhibit the COE level 
repulsion [10l|, and the distribution P(T) of the cavity 
transmission T at Ep reads 
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The thermopower Sk at a low temperature T (in units of 
TT^/cgT/Se, ks and e being the Boltzmann constant and 
the electron charge) is given by Sk = {dT/dE)/T where 
T and dT/dE are to be evaluated at E = Ep. These en- 
ergy derivatives can be obtained from the matrix Q, the 
distribution of its eigenvalues (the delay-times) taking a 
universal form for chaotic scattering: Assuming the WBL 
limit and a chaotic scatterer of size M ^ oo, the distri- 
bution of the inverse delay-times turns out to be given by 
the Laguerre ensemble from RMT [3, 0] • Introducing the 
two rescaled delay-times fj = tj/th (jh = l-Khj t^p be- 
ing the Heisenberg time with A p the mean level spacing 
of the cavity at E'f), tbe average eigenvalue density of 
the 2x2 matrix Q can be calculated using Refs. 
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while the distribution of the dimensionless thermopower 
CTfc = [l^p/2TT){dT/dE)/T can be found in Ref. in the 
form of a multiple integral: 

+ 1 oo oo +1 

PB{<7k) ^ j J J '^^^ J dT f{c,fi,f2,T) 



-10 



X 6 I CTfc -c(fi - f2)\l j; - 1 



(3) 




FIG. 1: (Color online) Chaotic scatterer of M = 7 sites (indi- 
cated in green) embedded in a Id tight-binding lattice (grey 
sites with nearest-neighbor hopping t — 1). The two end sites 
of the Id leads are shown in a darker grey. Upper right: Par- 
tition of the same infinite system used for deriving Eqs. ([5l-[6ll . 




Here /(c,fi,f2,T) = 

^(n,f2) 



,P(T)F(n,f2)F(n,T2) and 
cxp[-E-=i2?-]- The 
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charging effects inside the cavity are taken into ac- 
count 0, [3| in the Hartree approximation by the func- 
tion P(fi,f2): C being the capacitance of the cavity, 
F{fi,T2) ^ fi+T2 if eVC > Ap while F(n,f2) = 1 if 
eyC<^Ap. 

The distributions ([2]-|3l) were obtained in the WBL 
limit, an approximation where the real part of the lead 
self energy S(i?) and the energy dependence of its imag- 
inary part around Ep are neglected. This is justified in 
the bulk of a wide conduction band, but not near the 
band edges. Moreover, Hm was taken from the Gaussian 
orthogonal ensemble (GOE), a distribution giving rise 
to chaotic scattering only when M oo and Ap — >■ 0. 
These assumptions, leading to distributions can- 
not be made as £'_f approaches the edges of the conduc- 
tion band and the spectrum edges of the chaotic scat- 
terer. 

To study the delay-time and the thermopower distri- 
butions at the edges, we use a Id tight-binding lattice 
(hopping term t) embedding a scatterer of M sites, as 
shown in Fig. [T] S{Ep) is calculated from the scat- 
terer Hamiltonian Hm and from the exact expressions 
of the lead self-energies S(i?). The GOE distribution for 
S at Ep can be obtained by taking for Hm a Gauchy 
(Lorentzian) distribution of center Ep/2 and of width 
Tp = ty/l — {Ep /2t)2. The advantage of taking such a 



FIG. 2: (Color online) Average density of states per site (in 
unit of t~^), {p{E, Ep)) /M, as a function of the energy E (in 
unit of t), for 4 Cauchy ensembles (Hm £ C{M,Ef /2,Tf)), 
each of them giving rise to chaotic scattering at Ef ~ 
(brown line), Ef = 1 (red line), Ep = —1.5 (green line) and 
Ef = 1.9 (blue line). The circles indicate for each curve the 
value of the density of states at i? = Ep- The black dashed 
line is the semi-circle law (AfM)~^ (see Eq. (|12l) '). 



distribution for Hm is that it gives rise [17[ to chaotic 
scattering at Ep when M > 2, while one needs to take 
the limit M — >■ oo if Hm & GOE. A remarkable prop- 
erty of these Gauchy distributions is shown in Fig. [2] the 
average density of states per site at Ep varies on a semi- 
circle of center and diameter given by the center and the 
width of the conduction band, as one varies Ep. As Ep 
approaches the edges ±2i of the conduction band, Ep 
approaches the spectrum tails of the Gauchy scatterer. 
We will show that this gives rise to edge distributions 
for the time delay and for the thermopower which differ 
from the bulk distributions . The edge distributions 
turn out to be universal after an energy rescaling simi- 
lar to the one used by Tracy and Widom for the energy 
levels [i-Q. 

Scattering Matrix: Usually, the infinite system is di- 
vided into a scatterer and two attached leads, and the 
scattering matrix S is given [l8| in terms of the scatterer 
Hamiltonian Hm and of the lead self-energies S(-E). For 
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the Id-model sketched in Fig. [1] it is more convenient 
when M > 2 to divide the infinite system as indicated in 
the upper right corner of Fig. [T] a "system" made of the 
two sites (dark grey) located at the lead ends (Hamilto- 
nian Hq = V0I2), with the scatterer (Hamiltonian Hm, 
green) at its right side and the two leads (light grey) 
without their end sites at its left side. The scatterer and 
the leads are described by their self-energies [l^]- Using 
this partition, one obtains the matrix S at an energy E 
in terms on an effective 2x2 energy-dependent "Hamil- 
tonian" matrix H2{E): 



S{E) = -l2 + 2iT{E)A2{E) 
A2{E) ^ 



E12-Hq-H2{E)-11{E) 
H2{E) = W'^— — ^ „ W. 
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M 



(4) 
(5) 

(6) 



Here 1m is the M x M identity matrix and W is an 
Af X 2 matrix with Wi^j = 0, Wn = W22 = t- 
E = — 2tcosfc, Tj{E) ~ — te*'^l2 is the lead self-energy 
and r(£;) = isinfc = t^l- iE/2tf. Since Vq = 0, Ha 
disappears from Eq. (|6]). Using this expression for 5, 
dS/dE = 2id{T{E)A2{E))/dE and it becomes obvious 
that the distribution of dS/dE at Ep depends only on 
the distribution of the 2x2 effective Hamiltonian H2 (Ep) 
as Ep approaches the band edges ±2t (F^ — in this 
limit). We draw the conclusion that for having the dis- 
tribution of the delay-time and of the thermopower as 
Ep ±2t, one needs only to have the distribution of 
H2{Ep). Let us notice that when M = 2, it is simpler 
not to use the previous partition, and to put the 2x2 
Hamiltonian matrix H2 instead of H2 [E] in Eq. ([5]) . The 
matrices S given by the two methods are identical up to 
a phase factor e^"^^ . 

Chaotic Scattering and Cauchy Ensembles: An M x M 
Hamiltonian Hm has a Cauchy (or Lorcntzian) distribu- 
tion C{AI, e, F) of center e and width F if the probability 
of finding Hm inside the infinitesimal volume element 

1~T A'l 

fi{dHM) = Yli<j dHMaj is given by 

P{Hm) cx det {{Hm - el M 

For M = 2, it is straightforward to show that 

S{Ep) e COE^H2 e C{2,Ep/2,Tp). (8) 

Indeed, H2 and S are diagonalizable for M = 2 by the 
same rotation Rg and the eigenvalues of S are simply 
related to the eigenenergies Ej of H2: 
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e««. = -1 
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Ep/2- Ej +iTp 



; i = i,2. (9) 



Using Eq. ([T]) and Eq. (|9|) , one obtains the condition ([8|) . 

For M > 2, the condition 1^ for the effective Hamil- 
tonian H2{Ep) instead of H2 is sufficient and necessary 



for having chaotic scattering. Moreover, one can use two 
properties of the Cauchy distributions |17| . 

Hm e C(Af,e,F)^i/-/ e C{M,j^,%) 
Hm e CiM,e/D,T/D)^W^HMW eCi2,'^,'^) 



where D = + T^, to obtain the following result: 
Hm e C{M,Ep/2,rp) ^ S{Ep) e COE 



(10) 



since / D = 1 when D = Ep/A 

Chaotic Scattering and Semi-Circle Law: Hereafter, we 
study the distributions P{f) of the delay-time and P{ak) 
of the thermopower when iJjv/ G C{M,Ep/2,rp). In 
this case, S{Ep) € COE and the average energy level 
density of Hm at an energy E and for a Fermi energy 
Ep reads 



{p{E,Ep)) = - 



F, 



S{E-Ep/2)^+T% 



(11) 



where S — is the level spacing at E — 0. This gives 
for E = Ep 

Ap^ = {p{E = Ep,Ep)) = i^l - iEp/2t)^. (12) 

As shown in Fig. [21 the average level density per site at 
Ep of the scattering region must vary on a semi-circle 
as Ep varies inside the conduction band in order to give 
rise to chaotic scattering. This semi-circle is also the limit 
when M — >■ CX) of the average level density of a unique 
asymptotic Gaussian Ensemble where tr Hfj have a zero 
average and a variance = 26'^M/tt'^. Thus, chaotic 
scattering with {p{E = Ep,Ep)) becomes possible 
when Hm G C{M, Ep/2,rp) only if one of the energy 
distances = Ep ± 2t — 0. 

Energy Reseating: In our quest for universal asymp- 
totic distributions inside the conduction band (edges in- 
cluded), let us introduce the dimensionless parameter 
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The reasons for introducing a (i. e. for measuring the 
energy distance from the edges in unit of tM-^/"^) are 
twofold. 

Firstly, a appears as the relevant parameter in the 
study of the case M = 2, where the delay-time and ther- 
mopower distributions read [l9j : 



PM=2{T,a) 



4a 



(14) 



PM=2[o-k,a) = 2a\n j — j , (15) 

2Tra\ak\ 

where f = f — h/ {2Tpth) is a shifted delay-time. These 
distributions are valid in the whole conduction band 



for M = 2. Moreover, as pointed out after Eq. ([6]), 
they give the edge distributions for any Cauchy scatterer 
Hm e C{M,Epl2,Tp) of finite size M, in the hmit 
a ^ 0. 

Secondly, a is simply related near the band edges to the 
scale X leading to asymptotic universal level distributions 
at the spectrum edges of large Gaussian matrices. In 
particular, for an Hamiltonian Hm G GOE having the 
semi-circle density shown in Fig. [5] when M — >■ oo, the 
tail of its density {p'^^^{E)) /M takes the Tracy- Widom 
form outside the conduction band [20| 



1 



M 



exp[ 



r3/2l 



(16) 



in terms of the scaling variable x = e/{tM ^/'^). Since 
e=\E + 2t\ {e=\E- 2t\) ii E ^ -2t {E « 2t) respec- 
tively, one gets that x — (7ra)^/^ near the edges. It is 
clear from Eq. ((T5)) that the distributions corre- 
spond to the limit where x and hence a — > oo, since they 
were obtained in the asymptotic limit where p'^'~'^ {x) = 
outside the band. This limit is out of reach for M — 2 
(a < 2 /it). For finding the distributions ([2]-|3]) using our 
Id lattice model, we must consider larger scattering re- 
gions. This leads us to numerically check that our Id 
model does exhibit the bulk distributions ([2] - [3]) when 
a — >■ oo and the edge distributions (fT4l -[T5t when a — > 0. 
In addition, our results do confirm that the scaling lead- 
ing to Eq. (ITB)) leads also to asymptotic universal distri- 
butions for the delay-time and for the thermopower at 
intermediate values of a. 

Universal a-dependent asymptotic distributions: We 
study how the thermopower distribution P{ak,M,a) 
varies as one increases the size M of an Hamiltonian 
Hm € C(M, Ep/2,rp), taking Ep closer to the nearest 
edge for keeping a constant. To measure the size depen- 
dence of P{(yk, M, a), we introduce the parameter 



77(a, M) = 



/ dak\P{'Jk,M, a)~PB{<Jk) 
J dak\PM=2{a-k,a) - PsicFk) 



(17) 



where Psiok) [PM=2{(^k.a)] is given by Eq. © [Eq. ^] 
respectively. As shown in the inset of Fig. |31 P{ak , M, a) 
does reach an asymptotic limit when M ^ oo if a keeps 
a constant value. If this value is very small, 77 w 1 inde- 
pendently of M, confirming that Pa/=2 (c/c , a) gives the 
asymptotic edge distribution when a — > 0. For interme- 
diate values of a, the finite size effects on P{a-k,M, a) are 
noticeable when M is small, but become rapidly negligi- 
ble at larger M. We show also in Fig [3] how 77 depends on 
a as M increases. The curve obtained taking M = 100 
gives an excellent approximation of the asymptotic limit, 
and exhibits a crossover around a ~ 0.25 betwen the edge 
limit a — ?> and the bulk limit a — > 00. The same con- 
clusion can be drawn from the study of P{t, M,a). In 
Fig. m these asymptotic distributions (calculated taking 




FIG. 3: (Color online) 77 parameter as a function of a for 
M = 2 (dashed line), M = 3 (squares), M — 5 (diamonds) 
and M = 100 (circles). Inset : as a function of the size 
M of the cavity for a = 0.0025 (circles), a = 0.025 (down 
triangles), a = 0.125 (diamonds), a = 0.25 (up triangles), 
a = 0.5 (star) and a = 2.5 (squares). In both panels, full 
lines are guides to the eye. 



M — 200) are given for the delay-time and the ther- 
mopower for different values of a. The distributions 
P (fife, a), which are symmetrical around the origin, are 
shown for (jfc > only. The curves valid at the band 
edges and corresponding to the case M = 2 (Eq. (fH)) 
when a/Tp — and Eq. (fTS)) ') are also shown in Fig. |4l 
together with the bulk behaviors, given by Eq. ([2]) for the 
delay-time and by the numerical integration of Eq. (|3]) 
for the thermopower, taking F(Ti,f2) — 1. To show 
the convergence towards the edge limit a — > 0, we have 
multiplied the abcisses by a in the insets. One can see 
that the bulk distributions (|2l3p given in Refs. 0, 01 for 
chaotic cavities opened to leads via single mode quan- 
tum point contacts characterize also an infinite Id lat- 
tice embedding a Cauchy scatterer as a — )■ 00, while 
Eqs. (dH) and (fT51) give the edge distributions as a — 0. 
The change of the delay-time distribution as a function 
of a is striking. When a — >■ 00, the average delay-time 
(r) — irh/Ap and the fluctuations do not give rise to 
negative times. The scattering region is always attrac- 
tive. When a ^0, (r) ^ h/{2rp) (see Eq. ([HI)). This 
change of average delay time induced by a decrease of 
a is similar |21j] to the one induced in compound nuclei 
by a decrease of the density of resonance levels. Around 
this average value, there are huge symmetrical fluctua- 
tions which can yield positive or negative values for r, 
reflecting either an attractive or a repulsive character Q 
of the scattering region at the edges. 

In summary, we have obtained asymptotic edge dis- 
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FIG. 4: (Color online) Asymptotic delay-time (a) and ther- 
mopower (b) distributions for different values of a. Data are 
plotted as a function of f = t/th (a) and (b) in the 
main panels and as a function of a rescaled variable in the 
insets. Data for a = 2.5 x 10~^ (black dotted line in (a)), 
a = 2.5 X 10"^ (black dotted line in (b)), a = 2.5 x 10"^ 
(cyan line), a = 0.125 (magenta line), a = 0.25 (green dashed 
line) and a — 2.5 (blue dotted line) are shown. In the main 
panels, the data for a — 2.5 (blue dotted line) agree with 
the distributions (continuum blue line) expected in the limit 
a ^ cx) [Eqs. ((2]-E|)]. In the insets, the data for a = 2.5 x 10"^ 
(black dotted line in (a)) and a = 2.5 x 10"'^ (black dotted 
line in (b)) agree with the distributions (continuum red line) 
expected in the limit q — 5> [Eqs. (fTil - [T^ ]. 



tributions for the delay-time and for the low T limit of 
the thermopower using a Id model embedding a chaotic 
scatterer. Taking a Cauchy distribution for the scatterer 
Hamiltonian and the exact expressions for the lead self- 
energies, we have shown that chaotic scattering is asso- 
ciated to a semi-circle law for the average level density 
per site at Ep of the scattering region. This intrigu- 



ing coincidence between the asymptotic limit of a unique 
Gaussian ensemble and a continuum family of Cauchy 
ensembles has led us to use the same energy rescaling 
near the edges as the one leading to the Tracy- Widom 
expressions near the GOE spectrum edges. Doing so, we 
have shown that the delay-time and thermopower distri- 
butions converge towards universal distributions which 
depend on the scale a used for Tracy- Widom distribu- 
tions. In particular a — > oo corresponds to the bulk 
and a — > to the edges where we obtain analytic dis- 
tributions. We have numerically studied the universal 
crossover from the edges towards the bulk using the pa- 
rameter ri{a). 
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[1 
[2 

[s: 

[6: 
[7: 

[8 
[9 

[10 

[11 

[12 
[13 
[14 
[15 
[16 

[17 
[18 

[19 

[20 

[21 



S. F. Godjin, S. MoUer, L. W. Molenkamp, and S. A. van 
Langen, Phys. Rev. Lett. 82, 2227 (1999). 
S. A. van Langen, P. G. Silvestrov, and C. W. J. 
Beenakker, Supperlattices Microstruct. 23, 691 (1998). 
P. W. Brouwer, K. M. Frahm, and C. W. J. Beenakker, 
Phys. Rev. Lett. 78, 4737 (1997). 

P. W. Brouwer, K. M. Frahm, and C. W. J. Beenakker, 

Waves in Random Media (special issue on disordered 

electron systems 9, 91 (1999). 

F. T. Smith, Phys. Rev. 118, 349 (1960). 

C. A. Tracy and H. Widom, Commun. Math. Phys. 159, 

151 (1994)'. 

C. A. Tracy and H. Widom, Commun. Math. Phys. 177, 
727 (1996). 

P. J. Forrester, Nucl. Phys. B 402, 709 (1993). 

U. Doron, E ans Smilansky and A. Frenkel, Phys. Rev. 

Lett. 65 (3072). 

M. L. Mehta, Random Matrices 2nd edn (New York: Aca- 
demic, 1991). 

M.-J. Giannoni, A. Voros, and J. Zinn- Justin, Chaos and 
Quantum Physics (North-Holland, Amsterdam, 1991). 
R. A. Jalabert, J. L. Pichard, and C. W. J. Beenakker, 
Europhys. Lett. 27, 255 (1994). 

H. U. Baranger and P. A. Mello, Phys. Rev. Lett. 73, 
142 (1994). 

R. A. Jalabert and J. L. Pichard, J. Phys. I France 5, 
287 (1995). 

T. Baker, P. Forester, and P. Pearce, J. Phys. A 31, 6087 
(1998). 

P. W. Brouwer, S. A. van Langen, K. M. Frahm, 
M. Buttiker, and C. W. J. Beenakker, Phys. Rev. Lett. 
79, 913 (1997). 

P. W. Brouwer, Ph.D thesis, Leiden University (1997). 
S. Datta, Electronic Transport m Mesoscopic Systems 
(Cambridge University Press, 1995). 
A. Abbout, G. Fleury, J.-L. Pichard, and K. Muttahb, 
paper in preparation (2012). 

M. G. Vavilov, P. W. Brouwer, V. Ambegaokar, and 
C. W. J. Beenakker, Phys. Rev. Lett. 86, 874 (2001). 
V. L. Lyuboshitz, JETP Lett. 28, 30 (1978). 



